BOUNDS FOR THE LOSS PROBABILITY IN LARGE LOSS 
QUEUEING SYSTEMS 



VYACHESLAV M. ABRAMOV 

Abstract. Let G{si, 92) be the class of all probability distribution functions 
of positive random variables having the given first two moments gi and 92. 
Let Gi(x) and G2{x) be two probability distribution functions of this class 
satisfying the condition |Gi(a;) — G2(a:)| < e for some small positive value e and 
let Gi{s) and, respectively, G2{s) denote their Laplace-Stieltjes transforms. 
For real fi satisfying ^gi > 1 let us denote by jq-^ and -yQ^ the least positive 
roots of the equations z = Gi{fi — fiz) and z = G2(p — fJ-z) respectively. In 
the paper, the upper bound for {"fCi — IG2 I is derived. This upper bound is 
then used to find lower and upper bounds for the loss probabilities in different 
large loss queueing systems. 



1. Introduction 

In most of stochastic models studied analytically in the literature the probability 
distribution functions of their random characteristics are assumed to be known. In 
queueing problems, for example, the input characteristics are the distributions of 
interarrival and service times, and they are clearly described in the formulation of a 
problem. For example in the case of an M/G/1 queueing system, the arrival process 
is usually assumed to be Poisson of rate A, and service time distribution is assumed 
to be a given function B{x), with mean and other moments if required. This 
enables us to use the techniques of the Laplace-Stieltjes transforms or generating 
functions to obtain the desired output characteristics. 
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In practice, however, the distribution of an interarrival or service time is un- 
known. It can be only approximated by available information about that distribu- 
tion, and the accuracy of that approximation can be obtained by analysis of real 
observations. 

The problems of modeling, approximating and estimating the output character- 
istics of queueing systems are a well-known, and it is a well-established and distin- 
guished area of queueing theory. There is a wide literature related to this subject. 
To mention only a few papers that use different approaches, we refer Bareche and 
Aissani [Tl], Kalashnikov [12] and van Dijk and Miyazawa [20]. Bareche and Aissani 
[TT] used the strong stability method to study the error of approximation of GI/AI/1 
or M/GI /I queueing systems by that M/M/1, when the distribution of inter-arrival 
time or, respectively, service time is unknown but in the certain sense (that defined 
in that paper) close to the exponential distribution. Kalashnikov [12] studied sto- 
chastic sequences satisfying the recurrence relation Vn+i = -F(V„, where ^„ was 
a sequence of independent and identically distributed finite-dimensional random 
vectors. By replacing the original sequence {^n} by "perturbed" sequence {C^^}, 
under the assumption that a specially defined weighted distance between ^„ and ^'^ 
is given, that weighted distance between Vn and V^, where Vj[^i = has 
been studied. Van Dijk and Miyazawa [20j studied non-exponential queues such 
as GI/GI/l/n, M/GI/c/n and GI/M/c/n. For the GI/GI/l/n queueing system 
they demonstrated the influence of an error in the service time distribution on the 
resulting error in different performance measures such as the throughout of the sys- 
tem. They also established the error bounds for the throughout of the M/GI/c/n 
queue, and obtained similar results for the GI/M/c/n queue with a perturbation 
of the interarrival time distribution. 

In the present paper, we study the class ^(01,^2) of probability distribution 
functions of positive random variables having the given first two moments Qi and 
g2. We establish the bounds for the least positive root of the functional equa- 
tion z = G(/i — fiz), where G{s) is the Laplace- Stieltjes transform of an unknown 
probability distribution function G{x) belonging to the class G{gi,g2), and /i is a 
positive parameter satisfying the condition hqi > 1. The additional information 
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characterizing G{x) is that 

(1.1) 1C{G,F) —suplGix)- F{x)\ <€, 

x>0 

where F{x) is known probabihty distribution function of a positive random variable 
having the same two moments (i.e. belonging to the class G{qi,22) as well), and 
the least positive root, jp, of the functional equation z = F{^j, — fiz) is therefore 
known {F{s) denotes the Laplace-Stieltjes transform of F{x)). The metric /C(G, F) 
is known as the uniform (Kolmogorov's) metric (e.g. [T3], P^). 

The aforementioned bounds for the least positive root 7^ of the functional equa- 
tion z — G{fi — fiz) (or similar functional equations) are then used in asymptotic 
analysis of the loss probability in certain qucueing systems with the large number 
of waiting places. 

There are two areas of applications where these bounds are used. They are 
statistics of queueing systems and continuity of queueing systems. 

In statistical problems, the empirical probability distribution Gcmp(a;, A^) {N is 
the number of observations) is assumed to be known. If the number of obser- 
vations increases to infinity, then for any given positive value e the probability 
P |sup^>Q G'cmp(a;,^) — G{x) < e| approaches 1. 

More exact information about this probability is given by Kolmogorov's theorem 
(see KolmogorofF [15] or Takacs [19 , p. 170). Namely, 

lim p(.sup\G,mp{x,N)~G{x)\ < ^] ^ K{z), 

N~*oo {x>Q VN) 

where 

^ (-l)Je-2j'^', for z > 0, 
for z <0. 




So, the probability of 

(1.2) sup|G(a;)-Gcmp(a;,iV)| <e, 

x>0 

can be asymptotically evaluated when N is large and e is small. (For relevant 
studies associated with statistics (|1.2p or other related statistics see the book of 
Takacs 

The first and second moments of G(x) are usually unknown either. However, for 
large iV, they can be taken approximately to the empirical moments of Gomp(a;, N) 
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with some error. It will be shown below (see Theorem 12.31 rel. (I2.31|) and Remark 
12. 4p that if in relation (|1.2p the value e is small enough, then the bounds for the 
least positive root 70 are expressed via the first moment gi only. In this case 
the only error of the empirical mean is to be taken into account. Thus, in the 
motivation of assumption (|1.2p . the value e is assumed to be chosen such small that 
the probability of (|1.2p should be large on the one hand, and the bounds for the 
empirical mean should be small on the other hand. 

In continuity problems, we assume that the unknown probability distribution 
G{x) :~ P{( < x} with the expectation gi := j satisfies some specific properties 
such as 

sup \G{x)-P{C<x + y\C> y}\<e. 

x>0,y>0 

Then, according to the known characterization theorem of Azlarov and Volodin 
(see [TU] or [5]), we have 

sup|G(x) - (l-e-^^)l < 2e. 

a;>0 

For other related continuity problems see [6], where Kolmogorov's metric is used 
for continuity analysis of the M/M/l/n queueing system. 

The class of probability distributions functions G{0i,Q2) itself, i.e. without met- 
rical condition (|l.ip . has been studied by Vasilyev and Kozlov [21] and Rolski [17] . 
Rolski [T7] has established the bounds for the least positive root of the functional 
equation z = G(/i ~ /iz). 

In the present paper, we show that additional condition (jl.ip nontrivially im- 
proves the earlier bounds obtained by Rolski 17 . The new bounds have various 
applications. For example, the upper and lower asymptotic bounds can be obtain 
for the loss probabilities in M/GI/l/n, GI/M/l/n and GI/M/m/n queueing sys- 
tems with large capacity n as well as in many related models of telecommunication 
systems (see [5], [3], [1], [S], [7] and [S]). We demonstrate application of this theory 
for the GI/M/l/n queueing system with large buffer capacity n, for the M/GI/1 
buffer system with two types of losses [4 and then for the special buffers model with 
batch service and priorities J]. The last two of the mentioned applications have 
especial importance for telecommunication systems. We also establish new conti- 
nuity results for the loss probability in the M/M/l/n queueing systems with large 
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capacity n under special assumptions related to interarrival times. The continuity 
theorems for M/M/l/n queueing systems, where the buffer capacity n is fixed, 
have been established in [6]. Statistical analysis of M/GI/l/n and GI/M/l/n loss 
systems with fixed buffer capacity n based on Kolmogorov's statistics has been 
provided in [9]. 

The paper is structured as follows. In Section 2, properties of distributions 
belonging to the class G{0i, Q2) and satisfying additional condition (jl.ip are studied. 
Let Gi{x) and 6*2(2;) be arbitrary probability distribution functions of this class 
satisfying (|l.ip . i.e. Gi,G2 G ^(01,02), and /C(Gi,G2) < e. Denote by Gi(s) and, 
respectively, by G2(s) (s > 0) their Laplace-Stieltjes transforms. Let ■^Gi and 7G2 
be the corresponding solutions of the functional equations z — Gi{fi — fiz) and 
z — G2(/i — /iz) both belonging to the interval (0,1). (Recall that according to the 
well-known theorem of Takacs |18j . under the assumption /igi > 1 the least positive 
roots and 7G2 of the equations z — Gi{^ — ^z) and z — G2{n — ^z) are unique 
in the interval (0,1).) An upper bound for \jgi ~ 7G2I is obtained in Section [2l 
In Sections 3 and 4, applications of the results of Section [2] are given for different 
loss queueing systems. Specifically, in Section 13.11 lower and upper asymptotic 
bounds are established for loss probabilities in the GI/M/l/n queueing system as n 
increases to infinity; in Section [3?2l bounds for the loss probability in M/GI/1 buffer 
model with two types of losses, which has been studied in [1], are obtained, and in 
Section 13. 3[ bounds for the loss probabilities in the buffers model with priorities, 
which has been studied in [7], are established. In Section 21 the continuity analysis 
of the loss probability in the M/M/l/n queueing system is provided. The continuity 
analysis of Section 2] is based on the bounds obtained in Section [21 the results for 
the loss probabilities obtained in Section [3. II and characterization properties of the 
exponential distribution. 

2. Properties of probability distribution functions of the class G 

In this section we establish an inequality for \^Gi ^ IG2 1 for probability distribu- 
tion functions Gi(x) and G2{x) belonging to the class t/(0i,02) and satisfying the 
condition 



(2.1) 



sup|Gi(a;) - G2(x)| < e. 

x>0 
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We start from the known inequalities for probability distribution functions of 
the class ^(31,02)- Vasilyev and Kozlov ^ proved, 



(2.2) inf / e-'^'dGix) =c-'B\ s>0 

G6£;(fli,g2) Jo 

and 

(2.3) max / e-""dG(a;) = 1 - ^ + ^ exp -^s ], s > 0, 

GegiBuB2)Jo 02 02 V 0J 



where the maximum is obtained for 



(2.4) G{x) = G„,,^{x) - <^ 



0, if < < 0; 
1-^, ifO<t<fi^; 

82 ' — Bl 

1, ift>^. 

' — 81 

The lower and upper bounds given by (|2.2p and (|2.3p are tight. If 02 = 0? , then 
these bounds coincide. 

It is pointed out in Rolski T7] that (|2.2p and (|2.3p could be obtained immediately 
by the method of reduction to the Tchebycheff system if one takes into account 
that {l,i,t^} and {1, t, i^, e^**} form Tchebycheff systems on [0,oo). Rolski [17] 
has established as follows. For 70, the least positive root of the functional equation 
z — G{fi — fiz), it was shown 

(2.5) inf 7G = £, 

GGe(Bl,B2) 

and 

2 

(2.6) max 70 = 7G„„ = 1 + — (^ - 1), 

GGe(8l,B2) 02 

where £ in (|2.5p and (|2.6p is the least root of the equation: 

(2.7) z^e-mi+mi^^ 

The proof of (|2.5p and (|2.6p given in [T7] is based on the convexity of the function 
G{^ — /iz) — z. 

From p.2p and (|2.3p we also have as follows. Let Gi{x) and 6*2(2;) be arbi- 
trary probability distribution functions of the class G{3i,Q2), and let Gi(s) and, 

correspondingly, 6*2(5) be their Laplace-Stieltjes transforms (s > 0). Then, 

2 

(2.8) sup sup|Gi(s)-G2(s)| = 1-— . 
Gi,G2eg(8i,B2) s>o 32 
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Indeed, for the derivative of the difference between the right-hand side of (|2.3|) and 
that of (|2.2p we have 



This derivative is equal to zero for s = (minimum) and s = +00 (maximum). 

(The trivial case 02 — 2i, leading to the identity to zero of the right-hand side of 

(|2.9p for all s > 0, is not considered.) 

Therefore, from as well as from ([^ and we arrive at 

In turn, from (|2.5p and (j2.6p we have the following inequality for \jgi ^ 7G2I' 



The inequality (|2.10p follows from the results of Rolski Tf]. Under additional 
condition (|2.ip we will establish an improved inequality for \jgi ~ 702!- 

Prior studying the properties of the class of probability distribution functions 
G{Qi,Q2) under additional condition (|2.ip . note that inequalities (|2.2p . (|2.3p . (|2.5p . 
and ()2.6p hold true for a wider class of probability distribution functions than 
^(01 5 02)- We will prove that the above inequalities remain correct for the class 
of probability distribution functions m2)eM{gi 32) ^('^i' '^2)7 where the set of 
pairs {{mi, 1712)} contains the pair (31,32) (this set of pairs denoted by (01,02) 
will be defined below). 

Let m > i be such the boundary value, that the least root of the equation 

^ _ g-Aim-|-/jmz 

is equal to the right-hand side of (j2.6p . and let mi and 7712 are the values satisfying 

2 2 

the inequalities m < mi < gi, and ^ > ^ {ml < m2). Then, we have the same 
bounds (|2.2p and (|2.3p for the probability distribution functions and (|2.5p and (|2.6p 
for the roots 7^ but now for the wider class of probability distribution functions 
belonging to '?(0i,02) U Q{mi,m2). 

Indeed, for any mi satisfying the inequality m < mi < 0i, and any m2 for which 
— > — , according to [H] we have 

m2 — H2 ' ° ' ' 



(2.9) 





(2.10) 



|7G, -7G2l< 1 + 7^(^-1)- ^• 



02 



(2.11) 



inf / 

a(mi,m2) Jo 
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2 2 

and, taking into account that ^ > -Si and mi < ai together lead to — > — , we 
also have 

max / e dGix) = 1 ^ exp s 

G£gi7ni,m2) Jq 1712 "'-2 V "^1 

^ ' ? 9 / \ 

< 1 \ exp s , s > 0, 

02 02 \ 9i J 

where the equality in the right-hand side of the first line of (|2.12p is a replacement of 

the initial probability distribution function (given by p. 41) ) by another one, where 

the parameters 0i and 02 are correspondingly replaced with mi and m2. 

According to the result of Rolski [T7], we respectively have: 

(2.13) inf 7G = r>£, 

GeSimi,m2) 

and 

2 2 

(2.14) max 7^. = 1 + !Z^(r - 1) < 1 + ^(^ - 1). 

GeQ{mum2) 1712 02 

From the above inequalities of (|2.1ip - (|2.14p . one can conclude as follows. Let 

r Tfi^ 

A^(0i,02) = i ("^i'™2) : m < mi < 0i; — > — ; m^ < TO2 
L m2 02 

(Recall that m > ^ is such the boundary value that the least root of the equation 
z = Q-t^'^+t^-^^ is equal to the right-hand side of (|2.6p.) Denote 

Q{M)^ y e(mi,m2). 

(mi,m2)eA1(Bi,g2) 

Then we have the following elementary generalization of (|2.5|) and (|2.6p : 

(2.15) inf 7g = inf 7g — ^, 

GGe(gi,B2) Geg{M) 

and 

q2 

(2.16) max 7g = max 'fn^lG = 1 H ^ !)• 

Gee(gi,g2) Geg{M) 02 

2 

Notice that if toi = m, then we have £* = 1 + ■^^■(^ — 1), where i* is defined by 

2 

(|2.13p . On the other hand, according to (|2.14p we obtain ^ — 1, i.e. in this case 
TO2 = m^. Thus the set (01,02) and, consequently, the class G{M) are defined 
correctly. 
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We start now to work with ()2.ip . We have the following elementary property: 

/•oc />oo 

sup|Gi(s)-G2(s)| =:sup / c-'^'dGiix)- e-''^dG2(x) 

s>0 s>0 Jo Jo 

poo 

(o^7\ < sup / se^'^'^sup |Gi(y) - G2(y)| dx 

\^-^') s>OJo y>0 

v ' 

<e according to ( 12. H 

< £• 

Thus under the assumption of (|2.ip . the difference in absolute value between the 
Laplace-Stieltjes transforms Gi(s) and G2(s) is not greater than e. 
It follows from ([^17)) that 

(2.18) sup sup|Gi(s) - G2(s)| ei < e. 

Gi,G2ee(Bi,B2) s>0 
/C(Gi,G2)<c 

(We do not know whether or not the value ei can be found. However, the exact 
value of ei is not important for our further considerations. Relation (j2.18p will be 
used later in this section.) 

On the other hand, according to (|2.8p for two arbitrary probability distribution 
functions of the class G{.M) the difference in absolute value between their Laplace- 

2 2 

Stielties transforms is not greater than 1 — — . Therefore, if e > 1 — — , then 
the condition (jl.ip is not meaningful. Therefore, it will be assumed in the further 

2 

consideration that e < 1 — 

fl2 

The lemma below is the statement on the dense of the class GiQi,Q2)- 

Lemma 2.1. For any probability distribution function G{x) £ G{Qi,Q2) (0i 7^ 
Q2) there exists another probability distribution Junction G{x) £ ^^(01,52); which 
distinguishes from G{x) at least in one point, such that for any 5 > Q, 

/C(G,G) < 5. 

Proof. Under the assumption that the class G{qi, 02) is not trivial, i.e. 92 7^ , one 
can take two distinct probability distribution functions G{x) and F{x) of this class. 
For any pG (0,1), let Gp{x) ^ pF{x) + {\ - p)G{x). Apparently, Gp{x) e ^^(01,02) 
as well. Therefore, choosing p < |, by the triangle inequality we obtain: 

\Gp{x) - Gix)\ = \pF{x) + (1 - p)G{x) - G{x)\ < pF{x) +pG{x) <^-+^-=5. 

Hence, for p < | one can set G{x) — Gp{x). □ 
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An extended version of Lemma l2.1l is given in the following lemma. 

Lemma 2.2. Let (7711,7712) G M and {m[,m2) G M (7712 7^ ml, mi^ ^ {m[)'^), and 
let G{x) G ^7(7771, 7772). Then for any 5 > Q there exists a probability distribution 
function G{x) G (777 '1, 7773 ) such that sup^^g |G(a:;) — G{x)\ < 6. 

Proof. Assume that m'2 > 7^2- Then take a probability distribution function F(x) 
satisfying the properties: 

2777'i - (2 - S)mi 



27712 - (2 - <5)m2 



(2.19) / xdF{x) 

J —00 

and 

/oo 
x^dF{x) 
-00 



Then, the probabihty distribution function 

(2.21) G(x)= (^l-0G(x) + ^F(x) 

belongs to the class t/(mi, 7773), and, according to the triangle inequality 

\G{x)-Gix)\<^-G{x) + ^-F{x)<S. 

In the opposite case where m'2 < 7712 take the probability distribution function F{x) 
satisfying the properties 

, , , 27771 - (2 - S)m\ 
xdF{x) = \ 



2 , 2TO2 — (2 — (5)7772 



and 

x^dF{x) 

-00 

Then, instead of (|2.2ip we set 



G(2;) = (1 - G(2;) + ^-F{x). 

Apparently, G{x) G Q{m[,m'2), and, assuming that < F{x) < 1 for all a; > 
according to the triangle inequality we obtain: 

\Gix)-Gix)\<pix) + ^-F{x)<5. 

The lemma is proved. □ 
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With the aid of Lemmas 12.11 and 12.21 we wiU solve the following problem. Let 
Gi{x) and 6*2(2;) be two probability distributions belonging to the class Gi-M)- 
Under the assumption that sup^^Q 16*1(2;) — G2{x)\ < e we will find an estimate 
for the supremum of \"fGi — 7G2 1 (the supremum between the corresponding least 
positive roots of the functional equations z — Gi(/i — /iz) and z — 62 (/i ~ M^)-) 
Solution of this problem, in particular, addresses the case when the probability 
distribution functions 6*1(2;) and 62(2;) belong to the class ^(fli, 32)- In our analysis 
below the estimate of (|2.18p is used. 

The analysis uses Lemma [2.11 The Laplace-Stieltjes transform Gi(s) = e~^^^ 
contains only the parameter gi and does not contain the second one 02- Hence 
similarly to (|2.2p one can write 



(2.22) inf / e-"^dG(2;) = inf / e-'''"dG(2;) = e-^^^S s > 0, 
Gee(Bi,m2)Jo Geg{M)Jo 

where m2 is a fictive parameter, which is assumed to be unknown. Let us find 

this unknown parameter m2 in the Laplace-Stieltjes transform 62(5) = 1 — ^ + 

^ exp (-^s) taking into account that (cf. ^J^) 

(2.23) sup sup|Gi(s)-G2(s)| =ei. 

G2Ge(Bi,m2) s>0 
K:(Gi,G2)<c 

Relation (|2.23p holds true, because Gi(x) G G{.M), and according to Lemma [2.21 
for any e > there exists a probability distribution function Gi(x) e G{.M) such 
that |Gi(x) — 61(2;)! < e. On the other hand, the class of probability distribution 
functions ^(0i,TO2) is dense, so 61(2;) can be chosen belonging to the same class 
G {31,1^2) as the probability distribution function Gi(x). 

The real distance between the Laplace-Stieltjes transforms Gi(s) and 62(5) 
(Gi e G{gi,m2),G2 e G{Bi,m2)) is 

2 

(2.24) sup sup|Gi(s)-G2(s)| = 1- 

Gi,G2Ge(Bi,m2) s>0 ^2 

(cf. relation (12.8^ ). Therefore, equating the right hand side of (|2.24p to ei we have 

i ei. 



TO2 



and hence 



(2.25) m2 = -- 

1 — ei 
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The meaning of the parameter m2 given by ()2.25|) is as foUows. If the distance 
between two Laplace-Stieltjes transforms is ei in the sense of relation (|2.24p . then 
it remains the same for all distributions Gi{x) and G2{x) belonging to the family 
^^(flij 52), ^2 < 52 < £l2, where the class G{qi, ^2) is a marginal class of this family. 
In this case p.23|) can be simplified as 

(2.26) sup sup |Gi(s) — G2(s)| = sup sup |Gi(s) — G'2(s)| = ei. 

G2ee(gi,m2) s>0 G2ee(Bi,m2) s>0 

K{GuG2)<e 

Hence, in this case we have the bounds coinciding with the class of all distribu- 
tions of positive random variables having the moments mi = Qi and m2 — yz^, 
i.e. with the class Q ^gi, jz^^ ■ We also have as follows: 

sup |7Gi - 7G2I = 1 + — 1) --^ 



Gi,G2ee(^fli,j4^J i-Ei 
(2.27) =l + (l-ei)(^-l)-£ 

= ei - £ii 
<e-ee. 

Let us consider another case, where mi = gi — S > m, 6 > 0. Let Gi{s) = 

' ^\ / m2 m2 '^^ \ Si— 5 J 

parameter m2. In this case, 



with an unknown 



sup sup 101(5) - G2(s)| 

G2ee(Bi-<5,™2) s>0 
/C(Gi,G2)<c 

cannot be greater than e (see relations (|2.17p and (|2.18p ). 

For example, taking mi = m we arrive at Gi(s) = G2(s), and therefore 

sup sup |Gi(s) — G2(s)| = 0. 

G2ea(m,m=^) s>0 
>C(Gi,G2)<e 

For an arbitrary choice ofmi=gi — (5>m, one have 

sup sup |Gi(s) - G2(s)| = £2 < e- 

G2 65(01 -<5.m2) s>0 
'C(Gi,G2)<e 

(The exact value of £2 is not important.) In this case, similarly to (j2.25p 
(2.28) m2 = ii^, 

i — £2 
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and similarly to (|2^ . 

/o 9Q\ sup |7Gi - 7G2 I = £2 - £2^*, 

where i* is the solution of the equation z ~ i:^-i^ii:\-^)+p-{^i-^)^ _ \x readily seen 
that P > i. (The presence of positive (5 yields the value of the root of functional 
equation greater compared to the case where 6 is not presented (i.e. (5 0).) 
Keeping in mind that f > I and £2 < e, from (|2.29[) we have: 

Hence, from relations (|2.27p and (|2.30|) we arrive at the following theorem. 

Theorem 2.3. For any probability distribution functions Gi{x) and G2{x) belong- 
ing to the class ^(01,02) Bind satisfying condition (|2.ip we have as follows. 
//e < 1 - -Si, then 

■' fl2 ' 

(2.31) sup |7Gi -7G, I < e- e^- 

Gi,G2Ge(Bi,B2) 
/C(Gi,G2)<e 



Othe 



sup |7Gi - 7G2I = 1 + — (^ - 1) 
Gi,G2Ge(Bi,B2) 02 



'C(Gi.G2)<e 

where t is the least root of the equation 

2 = g-MBl+MBlz 

Remark 2.4. Theorem l2.3l is formulated for the class of distributions G{qi,Q2)- As 
it was mentioned, the parameters Qi and 02 are usually unknown. For practical 
applications, the errors for these parameters should be taken into account. Let 
us assume that ranges of these parameters are known. For example, (jri°'^°'' < 
01 < (^"PP°'' and 52°*°' < 02 < 52^'^°'^ are assumed to be satisfied with a given 
confidence probability P. It follows from Theorem 1 2 . 31 that if 52°™' > (5i^^°')^ and 
£ < I — ^^giowc? 1 then the lower bound i (the least root of equation (|2.7p ) for the 
least positive root 7^, should be replaced by the smaller value given by the least 
root of the equation 

This new value should replace £ in (|2.3ip to be used in real applications. 
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For a nontrivial class ^(01,02), and large enough volume of observations N, the 
above condition 52°*°' > {9i^^^'^)^ is natural. 

3. Asymptotic bounds for characteristics in large loss queueing 

SYSTEMS 

3.1. Loss probability in the GI/M/l/n queueing system. In this section we 
apply the results of Section[2]to large loss GI/M/l/n queueing systems (the param- 
eter n is assumed to be large). The results of this section are elementary. However, 
they serve as a basis for the analysis of the more realistic queueing systems, which 
are studied in Sections 13.21 and 13.31 The bounds for the loss probability obtained 
for this elementary system are then also used for a more delicate continuity analysis 
of the loss probability in M/M/l/n queueing systems in Section [4l 

Recall the known asymptotic result for the loss probability in the GI/M/l/n 
queueing system as n ^ cxi. 

Let A(s) denote the Laplace-Stieltjes transform of the interarrival time proba- 
bility distribution function A(x), let fi denote the reciprocal of the expected service 

time, let p denote the load, p = ~, , , which is assumed to be less than 1, and 

let a denote the positive least root of the functional equation z — A{p, — fiz). It has 
been shown in [3J that, as n ^ 00, the loss probability Pioss{n) is asymptotically 
represented as follows: 

(3.1) P\oss{n) = ^- +o(a ). 

I- p- p[l + p.A'{fi - fia)]a'' 

Notice, that the function ^'(z) = A{fj, — p,z) — z is a convex function in variable 
z. There are two roots z — a and z = 1 in the interval [0,1], and ^''(a) = 
~fj,A'{p ~ fia) — 1 > — 1. Therefore, according to convexity we have the inequality: 

(3.2) ^'{a) < 

a 

From p.2p we obtain: 

l + ^^A'{^^-^^a) > 

a 

and therefore 

(3.3) <i + ^,A'{n- ^ia) <1. 

a 
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Assume that A{x) G ^(01,02) is unknown, but the first two moments gi and 
02 are given. In this and the foUowing examples we do not discuss the statistical 
bounds for these moments such as those considered in Remark 12.41 So, all our 
examples are built on the basis of the moments gi and 02 only. 

Assume that Acmpix) is an empirical probability distribution function of this 
class, its Laplace-Stieltjes transform is Acmp(s), the root of the corresponding func- 
tional equation z — Acmpil^ — l^z) is a*, and according to available information, 
Kolmogorov's distance between Acmp{x) and A{x) is IC{A, Acmp) < e. 

Consider the case e < 1 — ^ (fli 7^ 02)- Since A{x) is unknown, A{s) will 
be replaced by Acmp(s) in p.3p . The numerator of the left-hand side of (|3.3p is 
replaced by the extremal element e^^^i^ which is not greater than that original. 
The corresponding denominator is replaced by (a* + e — e£), which is not smaller 
than that original a*. Assume that e is such small that a* — e + e£ > I and 
a* e - e£ < l-t- - 1). Then we have: 



(3.4) a*+e-e£ < ^ + M^Lp(M - < 1- 



Note, that the assumption on e under which (|3.4p is satisfied can be written as 
,^ ^ . /, 0? 02(1 ~a*)^ 0^(1-^) 



Using (|3.4[) . in the case of small e satisfying (|3.5p . according to Theorem 12.31 for 
n large enough we have the following two inequalities for lower P_(n) and upper 
P(ji) levels of the loss probability: 

(3.6) P(n) ^ (l-p)e--(a*-. + e£)'^ 



(3.7) P{n) 



(1 - p){q* d) - pe-i^si (^a* -e + e£)" ' 

(l-p)(a* + e-e£)" 



1 - p- p{a* + e - ££)«■ 

Therefore, for large n we have the following asymptotic bounds for Pioss{n): 

(l-p)e-^Bi(a* _e-he£)" 
(1 - p){a* +e- e£) - pe"™! {a* -e + tlY 

(3.8) < Pioss(n) 

^ (l-p)(a*+£-e^)" 
~ 1 - p - p(a* -I- e - e£)" ' 
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If e > 1 - then the terms (a* + e - ei) in ^1^) and should be 



replaced by these 1 + - 1) , and the terms (a* - e + ei) in dm and dSJ]) 
should be replaced by £. 



3.2. Losses from the M/GI/1 buflfer model. In this section we obtain lower 
and upper bounds for the loss probability of the following M/GI/1 buffer model [3]. 
Assume that messages (units) arrive in the buffer of large capacity N . Units arrive 
by batches, the sizes of which are independent and identically distributed positive 
integer random variables Vi with expectation c. In addition, the random variables ui 
are assumed to be bounded, i.e. p{;yi°™i' < Vi < = 1. Interarrival times of 

batches are exponentially distributed with parameter A, and the service (processing) 
times of these batches are independent and identically distributed random variables 
with the probability distribution function B{x) and expectation b. If upon arrival 
of a batch the number of units in the system exceeds the buffer capacity, then 
the entire batch loses from the system. In addition, there is probability p that an 
arrival batch of units does not join the system due to transmission error. In all other 
situations an arrival batch of units joins the system and waits for its processing. 

Assume that p = Xb > 1. Then, for the lost probability the following represen- 
tation has been derived in [4] (see relation (5.3) on page 757): 

(3.9) = . (P-l)+Ml + Ag(A-A/^)]E/^C(.) ^ 

P (p - 1) + [1 + AS(A - A/3)]E/3CW 

where B{s) denotes the Laplace-Stieltjes transform of the probability distribution 
function B{x), [3 is the least in absolute value root of the functional equation 
z = B{\ — Az), and 

C(iV) = sup |m : ^ < ivj . 

Notice that since p{f/i°wcr < < jyuppci j — 2, then the similar property for the 
random variable C,{N) is satisfied: P{C'°""(Af) < C(^) < C"pp"(^)} = 1, and, in 
addition, since as iV — > 00 

lim ^ = 1 

Af^oo N C 
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then, as -/V ^ oo, 

(3.10) E/?^(^) =/3^[l + o(l)]. 

Substituting ((3?T0)) into we obtain: 

P + (p-l)+p[l + Ag(A-A/?)]/3f , « 

(3.11) ttat = • — + o(p ■= ). 

P (p-l) + [l + AB(A-A/3)]/3f 

Similarly to (|3.3[) for the term 1 + XB'{X ~ A/3) we have the inequalities: 

(3.12) < 1 + Ai?'(A - A/3) < 1. 

Assume now that B{x) e ^2(51,52) is unknown, but with the first two moments gi 
and 02 are given, assume that -Bcmp(a;) is an empirical probability distribution func- 
tion of this class, its Laplace-Stieltjes transform is -Bcmp(s), the least positive root 
of the corresponding functional equation z = i?omp(A — Az) is /3*, and assume that 
according to an available information the Kolmogorov distance between -Bcinp(2;) 
and B{x) is /C(-B, -Bcmp) < £• Similarly to (|3.5p assume that e satisfy the inequality 

(3.13) e<min-;i , -, ; — 

^ ' I fl2 1-^ 02(1-^) 

Then similarly to (|3.4|) we have 

(3.14) (3*\T-ei - ^ + ^^-p(^ ^ - 

Under same assumption (|3.13p , using (|3.14p for large N we arrive at the inequalities 
for lower, tt^, and upper, ttjv, levels of this loss probability: 

p + p-l [p ~ + e - e£) + pe-^8i (/3* - e + e£)^ 
~ P (p-l)(/3*+e-e£) + (/3*+e-e^)f+i ' 

_ ^ p + p-l (p- l)(/3* +e-e^)+p(/3* +e- 6^)^+1 



P (p- l)(/3* +e- e^) + e-^Bi(/3* - e + e^)^ 
Hence, we arrive at the following statement. 

Proposition 3.1. Under the above assumptions given in this section for the M/GI/1 
fer model with large parameter N , in the case 



e < min <i 1 - , n e\ 

02 l-C 02(1-^) 
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for the loss probabilities ttn we have the inequalities: 



P ' (p-l)(/3*+e-e^) + (/?*+e-e^)f+i 

^ P + p-l {p~l){l3* +e-€i)+ p{l3* +e-e£)^+^ 
P (p- l)(/3* +e-e£) + e-^Bi(/3*_e f 

3.3. The buffers system with priorities. In this section we study the following 
special model considered in [7] (see also [S]). This model describes processing mes- 
sages in priority queueing systems with large buffers, and the effective bandwidth 
problem. Another interpretation of this system is a transportation system, in which 
vehicles pass by with some time intervals to pick up C passengers, at the most, with 
accordance of their priority status (different types of passenger are supposed to be). 
The formal description of the problem, given in terms of the buffers system with 
priorities, is as follows. 

Suppose that arrival process of customers in the system is a renewal process 
A{t), with the expected value of a renewal period j. There are I types of cus- 
tomers, and there is the probability pj > that an arriving customer belongs to 
the type J ^ -P'""''' = 1^ ■ Therefore, the time intervals between arrivals of the 
type j customers are independent and identically distributed random variables with 
expectation jpjj- 

Assume that for i < j, the customers of the type i have higher priority than the 
customers of the type j, so customers of type I are those of the highest priority 
and customers of the type I have the lowest priority. Assume that customers leave 
the system by groups of C as follows. If the number of customers in the system is 
not greater than C, then all (remaining) customers leave the system. Otherwise, if 
the number of customers in the system exceeds the value C, then customers leave 
according to their priority: a higher priority customer has an advantage to leave 
earlier. For example if C = 5, / = 3, and immediately before departure moment 
there are three customers of type I, three customer of type 2 and one customers 
of type 3 (i.e. seven customers in total), then after the departure there will only 
remain one customer of type 2 and one customer of type 3 in the system. Times 
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between departures are assumed to be exponentially distributed with parameter /x. 
Assume that ^ < 1- 

The buffer capacities for the type j customers is denoted N^^K Assume that all 
of the capacities N^^\ j = 1,2, ... ,^ are large enough, i.e. they are assumed to 
increase to infinity according to the rule that roughly is explained as follows. For 
specific numbers < ai < q;2 <...<«;< 1, the meaning of which is explained 
later, it is assumed that for any j < k 

(3.15) af 

where Nj :— ^ N^'^\ j — 1,2, ... ,1, is the cumulative buffer content of customers 

i=l 

of the first j types. 

k 

Let Pk '■— P'""''* be the probability of arrival of a customer of one of the first k 
types {pi = 1). Then the times between arrivals of customers, who are related to 
one of the first k types, k — 1,2, ... ,1, are independent and identically distributed 
with expectation 

Since ^ < 1, then pk = ^ < 1 for aU /c = 1, 2, . . . , /. Let Ak{x) denote the 
probability distribution function of an interarrival time of the cumulative arrival 
process generated by customers of the first k types, and let Ak{s) (s > 0) denote 
the Laplace-Stieltjes transform of Ak{x). For the Laplace-Stieltjes transform Ak{s) 
we have: 



(3.16) 

= PkAi{s) 



1 ~ {1 - Pk)Ai{s) 
Let Qffe denote the least positive root of the functional equation 

(3.17) z ^ Akifi ~ fiz^) 



(There is a unique root of this functional equation in the interval (0,1), see [?]■) 
Since pi < p2 < . . . < pi, then we also have ai < ^2 < • • ■ < It is shown in 
[7] and [S] that under assumptions (|3.15p . the loss probability of type k customers 
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is given by the asymptotic formula 



(3.18) (1 - Pfe)(l + afc + a| + . . . + - Pk[l + C^iA'^di - fia^)W^ 

+ ola 



) 

(The assumption = o (^a^''^, j < k, given in (|3.15|1 actually means that the 
losses of higher priority customers occur much more rarely compared to the losses 
of lower priority customers.) 

Our task is to find lower and upper bounds for nk ■ Note that asymptotic relation 
p.lSp is similar to that p.ip of the stationary loss probability in the GI/M/l/n 
queueing system with large n. The functional equation (|3.17[) is a more general than 
that considered before in Sections [T] and [2] (that functional equation is a particular 
case when C — I). For this functional equation, the lower and upper bounds are as 
follows. Let 01 := j and 02 denote the first and, respectively, the second moments 
of the probability distribution function Ai{x). (In the sequel, the notation 0i is 
used instead of y.) Then, for k = 1,2, ... ,1, from the representation of (I3.16P one 
can obtain the first and second moments of the probability distribution function 
Akix): 

/ xdAk{x) = — , 
Jo Pk 

and, respectively, 

2(1 -Pk)0l +Pfc02 



X dAkix) — „ 
Pk 

Furthermore, let £ denote the least positive root of the functional equation 
(3.19) z = exp' 



Pk 

(We use the same notation £ as it was used for the root of the simpler functional 
equation in Sections [T] and [21 because the consideration of the more general func- 
tional equation p.l9p leads to an elementary extension of the result of Rolski [T7] 
and consequently to elementary extension of the results in Sections [T] and [2]) . Fol- 
lowing this, we have: 



(3.20) 



inf Qffc = inf ua^ = £, 
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and 

(3.21) siipoffc = sup = 1 + -T- r^^- - 1), 

where a^i^ is the notation for the root of the above functional equation associated 
with the probabihty distribution Ak{x). (Along with the earlier notation a^, this 
notation is required for our purposes because it is spoken about the upper and 
lower bounds associated with the class of probability distribution functions defined 
in (|3.20[) . (|3.2ip and the equations appearing later in this section.) 

2 

Let us assume now that e < 1 — Trn P^, . Then according to the modified 

version of Theorem 12.31 related to this case we have the following: 

(3.22) sup \aA'^ - a A- \<t-d, 
A^'A^egy—, J 

/C(A',,A'fc')<e 

where aA' and aA" are the versions of ak corresponding the probability distribution 
functions A'f^{x) and A'^{x) of the class G ^il^Eh^+PhSl^ . 
Similarly to inequality p.3p . we have: 

(3.23) M^<i + cMUa*-M«^)<1, 

where A'f,{-) in (|3.23p denotes the derivative of ^fc(-)- 

Assume now that Akix) (^^, Hil:;£i:i|i±£Mi^ jg unknown, but the first two 
moments ^ and ^(^~P'')^fli+P'°B'' ^j-g given. Assume that Acrap,kix) is the empirical 
probability distribution function corresponding the theoretical probability distri- 
bution function Ak{x), and the Laplace-Stieltjes transform of Acmp,kix) is denoted 
by Acuip,k{s), s > 0. Let a^. denote the least positive root of the functional equa- 
tion z — Acmp.kifJ- — jJ.z'^). Assume also that according to available information, 
Kolmogorov's distance between Acmp,k{x) and Ak{x) is K,{Acrap,k, Ak) < e, where 
similarly to (j3.5p e is assumed to satisfy the inequality 

,^^,^,1 0? "fc-^ [2(l-p,)flf +p,02](l-«fe)-0?(l-^) 



2(1 - Pk)Ql + PkQ2 ' 1 - r [2(1 - pk)Q\ + PkQ2] (1 - t) 

Then similarly to (13. 4p we have 



exp 



(3.24) _^^<i + CMLp.fc(M-MK)'')<l- 
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Therefore, taking into account ()3.23p and ()3.24p for sufficiently large Nk we arrive 
at the following lower (denoted by E.k{^k)) ^^'^ upper (denoted by Tfk{Nk)) values 
for probability ttu ■ 

(f-pfe)exp (-Ml) (al-e + ee)^^ 



KkiNk) 



(1 - pk) E (a* + 6 - e£y - exp (-Mi) (a* - e + dY 
{l-puYal+e-eir" 



(1 - pk) E ("fc - e + ^^Y - PkH + e - e^)'^'" 



i=0 

Hence, we arrive at the following statement. 

Proposition 3.2. Under the above assumptions given in this section, in the case 
where 



e < niin < 1 



2(1-Pfc)fl?+Pfe92' 1-^' [2{1 - pk)3l + PkQ2]{l - i) 

for the loss probabilities nk, k ^ 1,2, . . . ,1, we have: 

(I-Pfc)exp(-Mi) („*_, + e^)^. 



(1 - pk) Eial + e - e£y - pk cxp (-Mi) (a* - e + ei^^ 

< TTfe < 



(l-pfe)(a^ + e-e^)^'= 



(1 - pu) E (a^ - e + e^)* - Pfc(a^ + 6 - et)^" 

i=0 

4. Continuity of the loss probability in the M/M/l/n queueing 

SYSTEM 

The results of Scction[2]enable us to establish continuity of the M/M/l/n queue- 
ing system when n is large. The continuity of the M/M/l/n queueing system was 
studied in 6J. In contrast to 6 where by continuity of M/M/l/n queueing system 
it is meant the continuity of a M /Gl /1/n queueing system, which is close to the 
M/M/l/n queueing system, in the present paper by continuity of the M/M/l/n 
queueing system it is meant the continuity of a GI /M /1/n queueing system, which 
is close to that M/M/l/n queueing system. Then, in the case when parameter n 
is large, the analysis becomes much simpler compared to the case when n is not 
assumed to be large. (In H] Conditions (A) and (B) mentioned below are applied 
to the probability distribution function of a service time.) 



BOUNDS FOR THE LOSS PROBABILITY 23 

Our assumptions here are similar to those of [6j. Let A{x) denote probability 
distribution function of interarrival time, which slightly differs from the exponential 
distribution E\(x) = 1 — e~^^ as indicated in the cases below. 

• Condition (A). The probability distribution function A{x) has the representa- 
tion 

(4.1) A{x) = pF{x) + (1 - p)Ex{x), < p < 1, 

where F{x) = Pr{C < x} is a probability distribution function of a nonnegative 
random variable having the expectation j, and 

(4.2) sup \Fy{x) - Fix)\ < e, e > 0, 

x^y>0 

where Fy(x) = Pr{^ < x + y\C > y}- Relation (|4.2p says that the distance in 
Kolmogorov's metric between F{x) and E\{x), according to the characterization 
theorem of Azlarov and Volodin ,ljOj (see also [B]), is not greater than 2e. 

• Condition (B). Along with (|4.ip and (14. 2p it is assumed that F{x) belongs 
either to the class NBU or to the class NWU. 

Recall that a probability distribution function S(a;) of a nonnegative random 
variable is said to belong to the class NBU if for all a; > and y > we have 
S(a; + y) < S(a;)S(j/), where S(a;) = 1 — S(a;). If the opposite inequality holds, i.e. 
S(x + y) > E{x)E{y), then S(a;) is said to belong to the class NWU. 

Under both of these Conditions (A) and (B) we assume that E^^ < oo is given. 

Under Condition (A), we have 

sup \A{x) - Ex{x)\ = sup \pF{x) - (1 - p)Ex{x) - Ex{x)\ 

(4.3) 

= psup|i^(a;) - Ex{x)\ . 

x>0 

According to the aforementioned characterization theorem of Azlarov and Volodin, 



sup|F(2:) - Ex{x)\ < 2e. 

x>0 



Therefore, from (14.31) we obtain 



(4.4) sup |^(a;) - Ex{x)\ < 2pe. 

x>0 
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We also have: 

2(1 -p) 



nOO 

/ x^dA{x) = pEC^ 
Jo 



A2 ■ 

Apparently, E^^ > (EC)^ — j7- Denote E^^ + p-, assuming that cr^ > 

Thus, it is assumed that Fj(^ > 

Now one can apply the estimate given by (|2.31[) in Theorem 12.31 to obtain 
continuity bounds for the loss probability in the case of large n. In this estimate, 
i is the least positive root of the equation z — exp (~ ^ + a-^)- It is not difficult to 
check that the least positive root of this functional equation is p = ^, and, because 
of the assumption > the value p is within the bounds: 

(4.5) i<p<i + _k^^e.i). 

So, keeping in mind the relation p.Sp for the loss probability, where in the given 
case a* is replaced by p, we arrive at the following statement. 

Proposition 4.1. Under Condition (A) and under the assumption cr^ > -p- the 
following inequalities for the loss probability, as n oo, hold: 

(l-p)e~p(p-2pei + 2pei£)" 

(1 - p)(p + 2pei - 2pei£) - pe^p (p - 2pei + 2pei^)" 

< Pioss{n) 

(I- p){p + 2pt2-2pe2iT 



< 



where 



and 



l-p-p(p + 2pe2 - 2pe2tT ' 

ei = min{p- ^,2pe(l - I)) , 



£2 = min ^ 1 + ^^ij-^ (^ - 1) - p, 2pe(l - ^) 



Under Condition (B) we have (|4.3p . where under the additional assumption that 
F{x) belongs either to the class NBU or to the class NWU one should apply Lemma 
3.1 of [5] rather then the characterization theorem of Azlarov and Volodin [TD], [5]. 
In this case we have 



sup \F{x) - Ex{x)\ < e. 

x>Q 



Therefore, from (14.31) we obtain 



sup\A{x) - Ex{x)\ <pe. 

x>0 
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In this case we have the following statement. 

Proposition 4.2. Under Condition (B) and under the assumption ct^ > -p- the 
following inequalities for the loss probability, as n ^ oo, hold: 

(1 - p)e~Hp - pes + pea^)" 

(1 - p) (p + pes - pes^) - pe" " (p - pes + pes^)" 

< Ploss("-) 

^ (l-p)(p + pe4-pe4y 
~ 1 - p - p(p + pe4 - pe4£)" ' 

where 

eg = min{p-£,pe(l-^)}, 

and 

£4 = min |l + - 1) - P,Pe(l - ^)} • 
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